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We propose a hybrid architecture in which an on-chip high finesse superconducting cavity is coupled 
to the lateral motion and spin state of a single electron trapped on the surface of superfluid helium. We 
estimate the motional coherence times to exceed 15 fis, while energy will be coherently exchanged with 
the cavity photons in less than 10 ns for charge states and faster than 1/is for spin states, making the sys- 
tem attractive for quantum information processing and strong coupling cavity quantum electrodynamics 
experiments. The cavity is used for non-destructive readout and as a quantum bus mediating interactions 
between distant electrons or an electron and a superconducting qubit. 



PACS numbers: 

The field of experimental quantum information process- 
ing has made significant progress in recent years. Many 
different physical implementations are being actively ex- 
plored, including trapped ions [1,2], semiconductor quan- 
tum dots [3, 4], and superconducting qubits [5, 6]. In par- 
ticular, the strong coupling to microwave photons possi- 
ble in circuit quantum electrodynamics (QED) architec- 
tures [7] has sparked interest in creating hybrid quantum 
systems capable of combining the advantages of different 
qubit implementations. In these proposals, a superconduct- 
ing transmission line cavity acts as an interface between su- 
perconducting circuits and microscopic quantum systems, 
such as polar molecules [8, 9], electron spins [10, 11] or 
ultra-cold atoms [12], typically with smaller couplings but 
much better coherence than superconducting qubits. Sin- 
gle electrons trapped above the surface of superfluid he- 
lium [13] might play a unique role as they can indepen- 
dently form a strongly coupled cavity QED system or act 
in concert with superconducting qubits. 

Interest in electrons on helium is motivated in part by 
their exceptional properties, including the highest mea- 
sured electron mobility [14] and long predicted spin co- 
herence times [15]. For these reasons the system was 
used in one of the first quantum information processing 
proposals [16]. The initial proposal focused on the mo- 
tional states of a single trapped electron normal to the he- 
lium surface [17], which promise long coherence times but 
have transition frequencies in the inconvenient range of 
100 GHz. Further, the electrons were to be detected de- 
structively. More recently it has been proposed to use elec- 
tron spins [15], and the possibility of moving electrons at 
MHz rates was demonstrated [18], but it was not clear how 
to best read-out or couple such spin states. 

Here, we address these challenges using the circuit QED 
architecture [7], and show that both the electron's motion 
and spin can be used to reach the strong coupling limit of 
cavity QED, where the coupling between the electron and 
cavity is larger than their decoherence rates, allowing for 



a wide variety of quantum optics and quantum informa- 
tion experiments. The quantized in-plane motion, parallel 
to the helium surface, can be engineered to have transition 
frequencies of a few GHz and is readily coupled to an on- 
chip cavity for non-destructive readout analogous to that 
used for superconducting qubits [7] or electron cyclotron 
motion in g-2 experiments [19]. The cavity also mediates 
interactions between individually trapped electrons allow- 
ing for multi-qubit gates similar to those demonstrated in 
superconducting systems [6]. In addition, the spin-photon 
coupling would be significantly enhanced by a controllable 
spin-orbit coupling. 

The trapped electrons can be considered as quantum dots 
on helium operating in the single electron regime. These 
dots would be sufficiently small (sub-micron) that the lat- 
eral spatial confinement and potential depth will determine 
the orbital properties. A Jaynes-Cummings coupling be- 
tween in-plane states and out-of-plane states in such dots 
was proposed recently [20]. The feasibility of creating such 
nano-scale traps is buoyed by a recent experiment which 
has detected single electron tunneling events [21]. How- 
ever, so far there have been no observations of either in- 
tradot quantization or spin resonance on helium. 

It is instructive to compare electrons on helium with 
semiconductor quantum dots. In most traditional two 
dimensional electron gases (2DEG's) such as in GaAs, 
the electrons form a degenerate gas with small effec- 
tive masses, renormalized g-factors, and strong interac- 
tions with the lattice. In particular the strong piezoelec- 
tric coupling leads to short coherence times for the mo- 
tional states (~ 100 ps) [22]. For this reason spin is typ- 
ically used [3, 4], but its coherence time can be strongly 
affected by nuclear spins [23]. In contrast, electrons on he- 
lium form a 2DEG at the interface between vacuum and 
superfluid, retaining their bare mass and g-factor. With the 
techniques described here, single electron quantum dots on 
helium promise some advantages over traditional semicon- 
ducting dots. We predict the decay of the orbital states to be 
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FIG. 1: (color) Top view of electrostatic electron trap. The 
ground plane and cavity center pin are shown in blue, while the 
trap electrode is magenta. The configuration of center pin and 
ground plane provide two-dimensional confinement. A DC volt- 
age, Vc is provided via a wire insulated from the resonator. Ma- 
nipulation and readout is performed via an RF voltage applied to 
the input port of the resonator with the modified signal measured 
by a cryogenic amplifier at the output port. 

10 6 times slower than in GaAs. Further, superfiuid 4 He has 
no nuclear spins (10~ 6 natural abundance of 3 He), lead- 
ing to long predicted spin coherence times[15], which are 
primarily limited by current noise in the trap leads. Per- 
haps most importantly, electrons on helium is a fascinating 
system where coherent single particle motion has not been 
accessible until now. 

An electron near the surface of liquid helium experiences 
a potential due to the induced image charge of the form 
V = -A/z, with A = e 2 (e - l)/4(e + 1) and e « 
1.057. Together with the 1 eV barrier for penetration into 
the liquid, the image potential results in a hydrogen-like 
spectrum E n = —R/n 2 of motion normal to the surface, 
with effective Rydberg energy R ~ 8 K and Bohr radius 8 
nm [13]. At the working temperature of 50 mK the electron 
will be frozen into the ground out-of-plane state, and the 
helium will be a superfiuid with negligible vapor pressure. 

With the vertical motion eliminated, the electron's lateral 
motion within an electrostatic trap could be coupled to the 
electric field of a superconducting transmission line cavity. 
As shown in Fig. 2, the cavity center-pin and ground plane 
form a split-guard ring around a positively biased trap elec- 
trode. We approximate the trapping potential in each of 
the lateral dimensions as being nearly parabolic, with level 
spacing ~ Tiu x ^ y . We assume a single electron in a high- 
aspect ratio trap so that the x and y motional frequencies 
are distinct, with u x < oo y . The Hamiltonian of the elec- 
tron near the potential minimum can be approximated as 

f 2 1 x 4 
H e = tt- + ^m e uj 2 x x 2 + ha— (1) 
2m e 2 6a% 

Here, a x = {Ti/meUJx) 1 ^ 2 is the standard deviation of the 
motional ground state wavefunction and a is the anhar- 
monicity. Because the trap is small and the potential must 
flatten at the outer electrodes, a < 0. The n to n+ 1 transi- 
tion frequency is io x _ n ~ u x ^ + (n+l)a. The electron mo- 
tion can be treated as a qubit when \ a\ is larger than the de- 
coherence rates. The scaling of the system parameters with 
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FIG. 2: (color) Side view of trap electrodes with energy levels 
and wavefunctions of electron motional state. The electron is 
confined to the surface of the helium film of thickness, d. The 
trap electrode (magenta) is biased positive relative to the ground 
and center pin (blue) of the CPW to laterally confine the electron. 
These electrodes form a confining potential which is harmonic 
to first order, but which flattens over the outer electrodes, giv- 
ing it a small softening anharmonicity. A sample potential and 
nearly-harmonic wavefunctions are shown. The spatial extent of 
the electron zero-point motion a x is small compared to the char- 
acteristic size of the trap w. To define a spin quantization axis 
a magnetic field in the x-direction is applied. To couple the mo- 
tional and spin degrees of freedom current is sent through the 
center electrode creating a z-field gradient within the trap. 



geometry (see Fig. 2) can be estimated analytically by ap- 
proximating the trap potential as V t cos(2nx/W). In this 
case oj x = 2Ti{eV t /m e W 2 Y' 2 , a = (2ir/W) 2 h/8m e , 
and Vt ~ V e e~ 27Td / w . Therefore one can tune the mo- 
tional frequency by adjusting the bias voltage, determine 
the anharmonicity by the trap size (confinement effects), 
and trade-off sensitivity in bias voltage for sensitivity to 
trap height (generally d ~ W so as to avoid exponential 
sensitivity to film thickness). 

The microwave environment and the trapping potential 
are simulated using Sonnet®and Maxwell®, respectively 
and then Schrodinger's equation is numerically solved to 
find the resulting wavefunctions for the geometry shown in 
Figs. 1 & 2. Using physically reasonable trapping param- 
eters: helium depth d = 500 nm, trap size W = 500 nm, 
trapping voltage V c = 10 mV results in a trap depth 
eVt/h ~ 20 GHz ^> k B T deep enough to prevent ther- 
mal escape, and a transition frequency oj x /2it ~ 5 GHz 
convenient to microwave electronics. The cavity can be 
represented by the Hamiltonian H r = Tioj r {a)a + 1/2), 
with uj r /2ir ~ 5 GHz close to the desired motional fre- 
quency. The electron's motion within the trap is affected 
by and induces an electric field in the microwave cavity. If 
the level spacing fiu x is in resonance with the energy of a 
cavity photon fvuj r , the two systems can exchange energy 
at the vacuum Rabi frequency, 2g = y/2ea x E /H where 
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E ~ 2 V/m is the zero-point electric field in the cavity. 
This yields a Jaynes-Cummings Hamiltonian of the joint 
system H = H e + H r + hg(a^c + ac^), where c is the 
motional quanta annihilation operator. The electron mo- 
tional states can be manipulated quickly due to the large 
coupling strength g/2ir = 20 MHz, a consequence of the 
large electron dipole moment ea x /\^2 ~ 2 x 10 3 Debye, 
and without exciting transitions to higher lateral states due 
to the anharmonicity a/27r ~ —100 MHz. 

In addition to the motional degree of freedom, the elec- 
tron carries a spin degree of freedom. The bare coupling 
of cavity photons to the spin is many orders of magnitude 
weaker than to the charge, but can be enhanced via con- 
trolled spin-motion coupling. A different mechanism of 
enhancement for semiconducting double-dots was pointed 
out in [24]. A spin-quantization axis is established us- 
ing a magnetic field in the x direction (Fig. 2). The Lar- 
mor frequency per unit field is approximately uj^/2-kB = 
2fi B /h ~ 2.89MHz/G. Niobium cavities have been 
demonstrated to maintain Q > 20, 000 in parallel fields 
of up to 2 kG, allowing Larmor frequencies of up to w L ~ 
6 GHz. Both the cavity and motion have fiui S> k B T so 
that they relax to the ground state. 

We propose to create a non-uniform z-field component 
with a gradient along the vibrational axis, d x B z , by passing 
a current through the center electrode (in the y-direction 
see Fig. 2). This leads to a new term in the Hamiltonian, 
H s = —2[i B s z xdB z jdx. The resulting spin-orbit inter- 
action provides an enhanced cavity coupling g s , mediated 
through the motional state, where 

a dB z gV2 
9s Mb x %oj x (1 — uj 2 l /uj1) 

This allows manipulation and readout of individual spins, 
as well as the use of coupling techniques developed for su- 
perconducting qubits [6]. Further, the coupling is propor- 
tional to the applied current, allowing the spin-cavity to be 
tuned in-situ on nanosecond timescales. For a 1 mA current 
500 nm away a dB z /dx ~ 8 mG/nm field gradient can 
be created. If oj l <C u) x these parameters give g s ~ 8 kHz 
whereas if uj x — u L ~ 30 MHz then the coupling can be 
made large, g s ~ 0.5 MHz. 

The current also creates a second-order variation in the 
a>component of B, leading to a new term in the Hamil- 
tonian, H sh = —fi B x 2 d 2 B x s x . If the constant magnetic 
field is applied along the y-direction, this term will lead 
to sideband transitions simultaneously changing the orbital 
and spin states for drives at uj± = lo x ± w L . These transi- 
tions can be used to manipulate, cool, and detect the spin 
using its coupling to the lateral motion. With such cooling 
it might allow one to use smaller spin frequencies. 

It is also important to consider decoherence of the mo- 
tional and spin states. The two major sources of noise are 
electrical fluctuations in the leads and excitations in the liq- 
uid helium. Here we present a short summary of these de- 
coherence mechanisms (also see Fig. 3). A detailed expla- 
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FIG. 3: (color) Decoherence rates of motional states as a function 
of the trap frequency due to interactions with bias leads, ripplons, 
and phonons. Rates are computed using parameters specified in 
the text at T = 50 mk. Solid lines are decoherence rates due to 
energy relaxation (ri/2), while dashed lines are dephasing rates 
(r^). Single ripplon relaxation rate and phonon dephasing rates 
are < 1 Hz. Spin decoherence rates are discussed in the text. 



nation of these mechanisms is presented in the supplemen- 
tary materials [25]. 

A motionally excited electron can relax radiatively via 
spontaneous emission directly into free space, through the 
cavity, or the trap bias electrode. The electron radiates very 
little into free space, both because it is small (a x <C A), 
and because the microwave environment is carefully con- 
trolled. In a perfectly symmetric trap, radiation through 
the bias leads would be suppressed by a parity-selection 
rule. We conservatively assume that the electron is dis- 
placed from the trap center by ~ a x , which gives a relax- 
ation rate ~ 1.6 x 10 3 s _1 . Though this mechanism is not 
expected to be dominant, it could be easily reduced signif- 
icantly by engineering the impedance of the trap bias lead. 
In addition, slow fluctuations in the trap electrode voltage 
(V e ) can deform the potential, changing the motional fre- 
quency and resulting in dephasing. This can occur from 
drift in the voltage source, thermal Johnson voltage noise, 
or local "1/f" charge noise. Drift slow compared with the 
experiment time is easily compensated. The thermal noise 
at 50 mK is quite small with dephasing rate < 100 Hz. Any 
charge fluctuations in the bias leads should be screened by 
their large capacitance to ground but even conservatively 
assuming an anomolously small capacitance, we estimate 
a dephasing rate 8 x 10 3 s _1 , which would not be the dom- 
inant decoherence rate. Noise from the cavity and ground 
plane electrodes should have less effect due to the symme- 
try of the potential. 

In addition to decoherence through the electrodes, the 
electron can lose coherence to excitations in the helium. 
Two major types of excitations are relevant: capillary 
waves on the helium surface, known as ripplons, and 
phonons in the bulk. The electron is levitated above the 
surface at height r B ~ 8 nm, which greatly exceeds the 
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height of the surface fluctuations, and therefore coupling 
to ripplons is small. The characteristic electron speed a x u> x 
significantly exceeds the speed of sound v s in He and the 
characteristic group velocity of ripplons. As a result, the 
rate of direct emission is suppressed and decay into rip- 
plons is dominated by second order processes in which 
two ripplons of nearly opposite momentum simultaneously 
interact with the electron. The allowed phase volume is 
limited by the condition on the total ripplon momentum. 
Thus the corresponding decay rate is small, estimated to be 
< 10 3 s- 1 (see Fig. 3). 

The most important mechanism related to helium exci- 
tations is decay into phonons. The coupling to phonons 
is reminiscent of piezoelectric coupling in semiconduc- 
tors. An electron creates an electric field that causes 
helium polarization, which in turn affects the electron. 
Phonons modulate the helium density and thus the polar- 
ization, which changes the electron energy. However, in 
contrast to semiconductors, where the typical piezoelec- 
tric constant is e pz ~ 10 14 e/cm 2 [26], its analog in He 
is ~ e(e — l)/47rr| ~ 10 10 e/cm 2 . Therefore coupling 
to phonons is much weaker than in semiconductors. The 
corresponding decay rate is ~ 3 x 10 4 s _1 (see Fig. 3). 

Besides decay, coupling to helium excitations leads to 
fluctuations of the electron frequency and ultimately to de- 
phasing. The major contribution comes from two-ripplon 
processes, since ripplons are very soft excitations with 
comparatively large density of states at low energies, so 
that they are excited even for low temperatures. However, 
because of the weak coupling, the dephasing rate remains 
small, ~ 2 x H^s" 1 for T = 50 mK (see Fig. 3). It 
also decreases rapidly as the temperature is lowered. An- 
other mechanism of dephasing are slow drifts in the helium 
film thickness, which change the trap frequency through its 
dependence on the height, d, of the electron. Fortunately, 
the cavity forms a liquid He channel [27] in which the film 
height is stabilized by surface tension, rendering it much 
less susceptible to low frequency excitations. 

The electron spin promises much longer coherence 
times, and when uncoupled to the charge the lifetime is 
expected to exceed seconds [15]. When the spin is coupled 
to the motion, it will also inherit the orbital decoherence 
mechanisms with a matrix element oc /i B d x B z a x /hu} x . 
These mechanisms can be further diminished by turning 
off the gradient field or changing the spin-motion detun- 
ing, to reduce the coupling. In addition to decoherence felt 
through the spin-orbit coupling, the electron spin can be de- 
phased by fluctuating magnetic fields. These can arise from 
Johnson current noise in the leads which would lead to de- 
phasing rates less than 1 s -1 . It is also possible that the 
spin will be affected by "1/f" flux noise [28], often seen in 
SQUID experiments. The trap involves no loops or Joseph- 
son junctions, so it is difficult to predict to the extent of flux 
noise in this geometry, however even a worst case estimate 
still yields a dephasing rate of only 200 s^ 1 [25]. 



In summary we use circuit QED to propose solutions to 
many of the problems associated with electrons on helium, 
developing the ability to manipulate and detect both the 
electron's quantized motion and its spin. Further, this ar- 
chitecture couples electrons on helium to each other and 
to other quantum systems via single microwave photons, 
creating a scalable architecture for quantum computing. 
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INTRODUCTION 



zero pressing field E± is 



These supplementary materials estimate the rates of 
relaxation and dephasing of the quantized electron mo- 
tion in a small electrostatic trap. A variety of de- 
coherence mechanisms are explored including, electric 
and magnetic noise from the electrodes as well as emis- 
sion and scattering of ripplons and phonons. We use 
DPS where we refer to Dykman et al, PRB 67, 155402 
(2003) [1]. The notations below are from this paper. We 
refer to Eq. (X) in this paper as DPS (X). A discussion 
of the relaxation and dephasing of the electron spin can 
be found elsewhere [2]. 



DEFINITIONS AND USEFUL RELATIONS 

In addition to the constants defined in the main text 
and DPS there are some definitions and relations that will 
be useful throughout these materials. The motional state 
of the electron is represented as \n x ,n y ) where n,{ x ^ 
is the quantum number in the x, y direction. The cor- 
responding transition frequencies are u>i and the local- 
ization lengths are a t = (h/mui) 1 / 2 , with i = x,y. In 
this supplement the electron is assumed to remain in the 
ground state of motion normal to the helium surface. 

With these notations, we have 



z (l\p 2 z /2m\l) z 



R 



mid) = |(i,o| e ^ r |o,o)r 



R = 



2mr| 



h x 158 GHz 



r B = h 2 /Am « 7.64 nm, 
A 



(6hc - l)e 2 
4(e hc + 1) ' 



(4) 



(1) 



where |1) 2 is the ground state wave function of motion 
normal to the helium surface. 

It is also helpful to include several properties of helium. 
The dielectric constant of helium is ehc ~ 1.057. The 
density of helium is p = 0.145 g/cm 3 , with surface tension 
a = .378 dyne/cm, and a dispersion relation for capillary 
waves (ripplons) is u) q = (crq 3 /p) 1 / 2 . For phonons the 
dispersion relation is ujq — Qv Sl with v s w 2.4 x 10 4 cm/s 
the speed of sound in superfiuid helium. 

For obtaining numerical results the following discus- 
sion will assume that the electron vibrational frequency 
is lu x w 5 GHz, the zero-point motion of the electron is 
a x rts 4.8 nm, and g rts 25 MHz. 



LIFETIME 

The general prescription of this section will be to write 
down Fermi's golden rule with the appropriate matrix 
element for the specific process. The matrix element will 
then be evaluated and integrated over the appropriate 
density of states. 



<frfc(q) = |(l,0|e ic > r |l,0)-(0,0|e^ r |0,0)| 



= - A (q x a x Ye-^^V 2 . 



Respectively, 



f dqg^q) = — — 

{2u) z J Aira x a y 

J^y /^(q) = — — 



8ira x a y 



(2) 



(3) 



The kinetic energy normal to the helium surface for 



Photon emission 

A vibrating electron can emit a photon via electric- 
dipole radiation into the vacuum. If the electron were in 
vacuum it would radiate at a rate 



,(v) 



2e2 
3 he 



2ira x 
A 



(5) 



where A = 2ttc/lo x is the wavelength of the vibrational 
frequency. If this were the only relaxation process the 
excited state would last for longer than 100 s. 

Of course, the electron is embedded in the electromag- 
netic environment created by the microwave resonator 
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and bias wiring. For a simple two-level dipolc coupled to 
a single mode cavity with frequency uj c for nonrcsonant 
coupling, g <C |w x — w c |, the spontaneous emission rate is 



r (v) 

1 1 



2 

g k 



(u x - LU C ) 2 



(6) 



where k is the mode decay rate. In the resonant case and 
in the bad cavity limit, where the cavity decay is faster 
than the coupling and the detuning, k 3> g, \uj x ~ w c|, the 
decay rate is r„ = g 2 /n. In the strong coupling limit, 
the rate of emission is = k/2 on resonance. It is 
also possible to model the case of a dipole in multi-mode 
cavity circuit [3]. 

In addition to the emission through the cavity, the elec- 
tron can also decay through the trap bias lead. This 
decay rate can be found by considering the effect of the 
Nyquist noise, 5y e (w x ) sa 2huj x Re[Z(uj x )} of the bias elec- 
trode on the electron. If there are no other fields, the 
bias electrode only couples to x 2 , and would not cause 
relaxation. However any displacement Ax from the bias 
electrode null, due to stray or intential DC fields from 
other electrodes will add a x interaction and will open 
a channel of relaxation. In this case the coupling to the 
electrode potential fluctuations 5V e is given by hSV e , with 
h = exdE x /dV e , where E x is the electric field on the 
electron in the x direction due to the bias electrode. The 
decay rate is then 



F (el) _ 1 



(0,0|fc|l,0) S v ,(u x ). 



(7) 



For small Aa; compared to the inter-electrode distance, 
we have eE x = m e uj 2 Ax = hw x Ax/a 2 . Substituting 
these results into Eq. 7 and allowing for dE x /dV e — 
E x /V e , one obtains 



.(el) 



Rc[Z(u x )} 
h/e 2 



~eV, 



u x 



(8) 



Even assuming that the electron is displaced significantly, 
Ax « a x and Ke[Z(uj x )] ss 50 fl, this relaxation is quite 
small « 1.6 x 10 3 s _1 and should not be a limiting 
factor. If this becomes a hindrance in the future it can 
be reduced by several orders of magnitude by changing 
Rc[Z((x) :I ;)] using a resonant structure. 



Two-ripplon scattering 

Single ripplon scattering is exponentially suppressed 
by the mismatch of the size of the electron wave func- 
tion to the ripplon wavelength at the same energy. The 
dominant decay process is one in which two ripplons are 
emitted in nearly opposite directions, each with approx- 
imately half of the electron energy Huj x . 



The decay rate due to such a process can be estimated 
by applying Fermi's Golden rule to DPS(ll) 

r ( 2r) = 2n ^ | ( o,0|e qi U e 4(qi+q2) ^ qi , q2 |l,0) 



qi,q2 



xS(Hlj — huj qi — huj qi )(n 



+ 1), (9) 



where £ q is the helium surface displacement due to a 
ripplon with wave vector q and n q = n(cj q ), with n(oj) = 
[exp(fiu! I kftT) — 1] . This calculation can be simplified 
significantly by noticing that, because of the Gaussian 
form of the matrix element Eq. (1), |qi + q 2 | < l/a^- 
Since q ICS found from condition 2w groB = lo x is w 4.6 x 
10 6 cm and q Tes a x sa 30, we approximate |qi| w | q 2 1 



I : ■ strength of the direct two-ripplon coupling Vq iq2 is 
a sum of contributions from the kinetic, or inertial term 
(from electron accommodating to the curvature of the 
barrier on the helium surface), and a polarization term, 
which comes from the change of the image potential due 
to the curvature of the surface. In Eq. 9 and below it 
should be assumed unless otherwise stated that T^ iq2 is 
the matrix element of the interaction with ripplons or 
phonons on out-of plane wave functions |l) z . A simple 
calculation shows that the main contribution to scatter- 
ing comes from the kinetic term, see DPS (11)-(13), 



-i?qiq 2 . 



(10) 



A correction from the polarization term, DPS (13), is 
comparatively small. The corresponding matrix element 
of the interaction is 



Az~ 3 [2 - x 2 K 2 (x)] ss 4.6R/r% 



(11) 



with x = q les z. The ratio of this term to the kinetic- 
energy term is 40%. The remaining term in the polar- 
ization coupling, DPS (13), is even smaller and partly 
compensates the above term. The contribution of the 
polarization coupling decreases, relatively, when we ap- 
ply a pressing field. 

Using this approximation to evaluate the matrix ele- 
ment and converting the sum to an integral yields 



(2r) 



J w; 



16ttV 
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J dqgi(q) 



(q 2 R) S(huj - 2hw q ){2n q + 1) 



R 



n 2 / 3 , , 7/3 



(12) 



192™^ 2 1 /3 cr 8/3 
This rate is small, <~ 450 s _1 . 

Renormalization due to single-ripplon coupling 

The single-ripplon coupling leads to renormalization 
of the two-ripplon coupling. The renormalization leads 
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to the replacement of the matrix element of two-ripplon 
coupling 



V^(q) = (Mle^Vqk). 



(13) 



Here, \fi) are electron states enumerates both in-plane 
and out-of-plane states), are the state energies. We 
have disregarded the dynamics of ripplons, i.e., we as- 
sume that the change of the electron energy in a virtual 
transition is much larger than the ripplon energy change. 
The explicit form of the coupling matrix element (an op- 
erator with respect to motion normal to the surface) Vq 
is given by DPS (10). We need virtual transitions into 
such states /i that the matrix elements (//y | exp(zqr)|0, 0) 
is of order 1 for q <~ q rcs , where is the in-plane com- 
ponent of the wave function \fj). This condition leads 
to typical « h 2 q 2 cs /2m. Such energy is very large, 
~ 12R. The terms in Vq that are oc q 2 are also large for 
such q. They lead to renormalization [4] of the matrix 
elements of direct two-ripplon coupling, reducing them 
by a factor oc (gr B ) _1 which is ~ 0.3 for q = q rcs . The 
polarization coupling is also reduced. Therefore, the rate 
of two-ripplon decay is reduced down to < 10 3 s _1 . 



Phonon scattering 

Another way that the electron motion can relax is 
through phonon emission. The electron can launch a 
phonon by two mechanisms. In the first the coupling is 
mediated by modulation of the dielectric constant along 
with the density wave. In the second, the electron can 
couple to the surface displacement in much the same way 
as in the ripplon case but now launching most of the en- 
ergy normal to the surface. 



Phonon scattering: dielectric constant modulation 

For frequency oj x = 5 GHz and phonon sound velocity 
v s = 2.4 x 10 4 cm/s, we have the normal to the surface 
component of the wave vector of the resonant phonon 
{Qz)r = Wj/fs ~ 1-3 x 10 6 cm -1 , whereas the typical in- 
plane wave number q x = (mujx/h) 1 ^ 2 w 1.6 x 10 5 cm -1 . 
This shows that the phonons involved in inelastic scat- 
tering propagate almost normal to the surface. We will 
modify DPS (25) - (27) to write the interaction in terms 
of the velocity potential 

00, z ) — 0Q (cqQ, + cL qQ J e iqr sin(Q z z), 



Q 



4>a = (h Vs / P VQ) 



1/2 



Q = (q,Qz), 



(14) 



where c q q z is the annihilation operator of a phonon with 
the wave number (q, Q z ); r,q are two-dimensional vec- 
tors and Q z > 0. Using the Fermi golden rule, for 
n(co x ) <C 1 the decay rate can be written as 



r< d > = 



2tt 



Q 



V, 



(d) 



Q 

(d) 



xS(huj x — Hluq^), 

^,<i|t;W|i>„ 



(15) 



I 



OO 



dz'(z + z'y 1 sin(Q z z')^i[<z(z + z')]. 



Here, ujq h ' is the phonon frequency. 

Using the explicit form of the ground-state wave func- 
tion and taking into consideration that the typical in- 
plane wave numbers are small, q <C l/r\B, so that in 
Eq. (15) Ki[q(z + z')] w l/<z(z + z') (this is an overes- 
timate), one can write the expression for the decay rate 
due to modulation of the dielectric constant e as 



= 32R 2 uj x / (71 pi/gha x a y ) 
2 



dzi d,Z2 ( — 

\zi+ z 2 



e 2zi sm(Q z r B z 2 ) 



(16) 



The dimensional factor in front of the integral in 
Eq. (16) is 4.3 x 10 6 s" 1 . With (Q z ) r r^ 0) = 1.01 this 
gives the overall rate of <~ 2.7 x 10 4 s" 1 . 



Phonon scattering: surface displacement 

The contribution to the decay rate comes also from the 
displacement of the helium surface, which can be consid- 
ered a free surface for high-energy phonons propagating 
almost normal to it. This is much like the ripplon case, 
but here a single-phonon decay is allowed because the 
2D momentum is conserved while the excess energy is 
dumped into the normal component. The decay rate can 
be expressed as 



r i s) = ^EK°'°i^ e<q "^i 1 ' >l 



x5{fkj - ftw$ h) ) 



(17) 



Here, U q is the diagonal matrix element of the coupling 
operator U q , which is detailed in DPS (9) - (10), calcu- 
lated on the wave functions of the ground state of motion 
normal to the helium surface, 



Vq - ,(l|Vq|l) 2 - eE ± + Aq 2 (v pol ), 



(18) 



where E± is the electric field that presses the elec- 
tron against the surface and (v po \) = z (l\(qz)~ 2 {l — 
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qzK x (qz)\l) z . Then 



hpv s u x (2tt) 2 



The contribution from the pressing field is 



1 1 — 



2nhpi/ s uj x a x ay 



(19) 



(20) 



300 V/cm this 



For a comparatively strong field Ej_ 
rate is 8.6 x 10 3 s" 1 . 

We now consider the contribution to 
term (vpoi). For typical q < l/a x we have (i> po i) ~ 
— ln(grs)/2. Using this approximation, we obtain for 
the corresponding contribution 



r^ s) of the 



p(pol) _ 



A 2 



- / dqq 4 g l (q)\(iy pol )\ 1 

X J 



~ 48J?2r| „ [ln(r B /a x )/2] 2 . (21) 



Here, we assumed that a 2 



The last factor is some- 



what smaller than the numerical result for actual rB/a x . 
For a numerical estimate we set it equal to one. This 
gives 1.8 x 10 3 s _1 . There is an interference term of the 
two last expressions, but it is smaller than their sum. 

Therefore, polarization phonon scattering is the dom- 
inating mechanism of inelastic scattering due to excita- 
tions in liquid helium. 



low frequency noise on the bias electrode will result in 
electron phase noise. Noise on the bias lead can arise 
from slow drift of the voltage source, thermal fluctuations 
(Johnson noise), or anomalous local sources (1// noise). 
For typical electrode geometries the transition frequency 
scales like uj x cx V e 1//2 . Small changes in V e cause a fre- 
quency shift Su! — —U! x SV e /2V e . Typical long term drifts 
for precision voltage sources are 5V e /V e « 10~ 6 , or a drift 
of about 5lj x w 8 x 10 3 s _1 . Because this drift is typically 
on hour timescales it should be easily compensated by 
measuring the transition frequency and readjusting the 
voltage. For thermal noise [5] 



r («) = 



Svj2 



Lj 2 k B TRe [Z(w)] 



4K 2 



(22) 



If superconducting leads are used then there is no local 
Johnson noise but noise can still be coupled in from dis- 
sipative sources off chip. For the simulated bias voltage 
V e w 10 mV and conservatively assuming an environmen- 



tal impedance of Z — 50 f2, Eq. 22 gives T 



90s' 



Usually the dominant source of electrical dephasing is 
1/f noise which is thought to arise from mobile charges 
in the substrate. This type of noise spectrum does not 
lead to a simple dephasing rate. The diffusion of phase 
can be estimated as [5] 



Lifetime summary 

There are several contributions to the lifetime of the 
excited vibrational state of the electron. There are ef- 
fects which could modify these results (though the esti- 
mates here are conservative), but in the current analysis 
the dominant mechanism appears to be phonon emission, 
coupled mostly via modulation of the dielectric constant. 
It gives the lifetime of the vibrational state T\ « 35yus. 



DEPHASING 

In addition to loss of coherence of the electron due 
to decay processes, excitations in the environment at 
low frequencies can cause variations in the transition fre- 
quency, leading to dephasing of superposition states. In 
this section we treat dephasing from voltage noise, vari- 
ations in helium level, and state dependent quasi-elastic 
scattering of ripplons off the electron. 



Voltage noise 

The depth of the potential and thus the electron tran- 
sition frequency depends on the trap bias voltage. Any 



<[</>(*) -0(O)] 2 > 



3uj x 
dV e 

W 2 



5y e (lHz) ln(0.4/f m t)t 2 
S v .(lHz)In(0.4/f m t) (23) 



where /~ x is the total averaging time, and t is the pre- 
cession time of single measurement. The dephasing time 
for uncompensated 1/f noise can be estimated as the 
time for {\4>(t^) — (f)(0)] 2 ) = tt. Local charges should be 
screened by the bias electrode. To get an estimate of 
their effect we assume that Sy e = S q /C g and that the 
effective capacitance is only C e g w 1 fF (though it should 
be much more than that) and the "typical" charge noise 
of 5,(1 Hz) w lO^c/Hz 1 / 2 . With these assumptions 
t^ 1 = 8 x 10 3 s _1 . Even these conservative assumptions 
lead to rates smaller than the relaxation rates, which can 
be seen as a result of the electron being physically sepa- 
rated from the electrodes where the charges may fluctu- 
ate. 



Fluctuations in helium level 

The potential energy landscape is created by electro- 
static gates beneath the helium surface. Any fluctuations 
in the thickness of the helium above the electrodes will 
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change the effective voltage seen by the electron result- 
ing in changes to the transition frequency. While quan- 
tized excitations in the helium level will be treated in the 
next section, here we explore the effect and susceptibil- 
ity to slow changes in the level of liquid helium in the 
reservoir due to fluctuations in temperature, or external 
vibrations. This type of fluctuation is unlikely to cause 
dephasing during the lifetime of the qubit and should be 
susceptible to spin-echo techniques, but if unchecked can 
lead to slow drifts in the frequency, making it difficult to 
bias the electron reliably. 

For the geometry simulated in Fig. 2 the transition fre- 
quency has a sensitivity to small fluctuations in the mean 
helium film thickness of du> x /dd ss lOMHz/nm. Fortu- 
nately superfluid has two properties, the formation of a 
Van der Waals film, and ideal capillary action that help 
stabilize the film thickness. Exploiting these properties, 
the electron can be trapped above a channel [6] in which 
the thickness is determined by the geometry and surface 
tension. The trap is supplied by the Van der Waals creep 
film from a reservoir located well below the trap [7] , sig- 
nificantly reducing its sensitivity to fluctuations in the 
reservoir height. To accomplish this we design the guard 
ring of the trap to be much thicker than the bias elec- 
trode, forming a channel of the desired height. The su- 
perfluid obeys Jurin's law, R c = 2a/pgH causing it to 
fill the trap if it is smaller than R c , the capillary radius 
of curvature [6] . For H = 5 mm, the capillary length is 
R c s=s 53 /im and the trap is well filled, with sensitivity 
reduced by duR c dR c d = w/H w 10~ 4 . 

Direct two-ripplon coupling 

An important contribution to dephasing is due to ther- 
mally excited ripplons. For T= 100 mK the typical 
wave number of a thermal ripplon (hcj qT = k^T) is 
q T sa 4.1 x 10 6 cm" 1 . Therefore q T r B w 3.1. For T=50 
mK we have qx ~ 2.5 x 10 8 cm -1 and qrvs ~ 1.96). 
The dephasing rate is given by the difference of the diag- 
onal interaction matrix elements on the wave functions 
of the excited and ground vibrational states. The overall 
interaction matrix elements need to be projected on the 
ground state of out-of-plane motion, that is the relevant 
interaction can be written as [cf. DPS (32)] 

H t = E E ^•(qi,q2)^ 1 & q2 |j,0)0',0| (24) 

j=l,2qi,q 2 

with 

^(qi,q 2 ) = (h/ p)(q 1 q 2 /uj qi uj q2 ) 1/2 V cll -q 2 
x(j,0|exp[i(qi -q 2 )r]|j,0). 

A useful relation for the following calculation is Eq. (2). 
For qxTB S> 1 the major direct two-ripplon coupling 
is the kinematic coupling, with the projection given by 



Eq. (10). With account taken of DPS (32)-(33), the de- 
phasing rate given by this interaction has the form 

r (2r) _ 

^E ( t \ ) ( c ii f i2) 2 .gph(qi +q 2 ) 

xn qi (n q2 + 1) 6 (u qi - oj q2 ) (25) 

(S is the area). For q T ^> l/a x ^ y , we can calculate this 
expression using that |qi — q 2 1 "C gi, 2 ~ It- The calcula- 
tion is further simplified if we set a x = a y , in which case 
we can average over the directions of qi.2, which leads 
to the replacement g p h(q) — > (3/32)(ga x ) 4 exp(— q 2 a 2 /2). 
We can then integrate first over q 2 for given qi by writ- 
ing q 2 = g'qi + g"(qi x z) (hat indicates a unit vector). 
The inequality |q x — q 2 | <C gi ;2 shows that \q'\ >• \q"\. 
To leading order in (a^g^)" 1 we have uj q2 « co q i and 
qi q 2 ~ <z"(qi x z). Integration over q',q" then be- 
comes simple, and the remaining integration over qi is 
then reduced to integration of the weighted temperature- 
dependent factor n qi (n qi +1). The overall result is 

r (2r) _ (2^) 1/2 pR 2 . . .3 , . 

For 100 mK this gives - 1.4 x 10 4 s" 1 . 

For T~ 100 mK the polarization part of the direct 
two-ripplon coupling is small compared to the kinematic 
part. However, its role increases with decreasing tem- 
perature, since the characteristic ripplon wave number 
qr decreases. For qi « q 2 = q, the matrix clement of 
the polarization coupling — z (l|Az~ 3 [2 — (qz) 2 K2(qz)]\l) z 
changes from w — 0.9Rq 2 to w — 0.6Rq 2 for q increasing 
from 0.3/rs to 1.2/rg. It is thus comparable with the 
kinematic coupling Rq 2 and has opposite sign. There- 
fore the overall contribution to the dephasing rate from 
direct two-ripplon coupling is significantly less than from 
the kinematic coupling taken alone, which from Eq. (26) 
is - 1.7 x 10~ 3 s- 1 for T=50 mK. 



Renormalization due to single-ripplon coupling 

The kinematic and polarization single-ripplon cou- 
pling taken to the second order of the perturbation the- 
ory renormalizes the matrix elements of the direct two- 
ripplon coupling. For qi^B ^ 1 it reduces these matrix 
elements by a factor ~ l/qipfB [4]. This leads to a sig- 
nificant reduction of the dephasing rate for T=100 mK, 
since q^s ~ 3 for such temperature. On the other 
hand, for lower temperatures down to T=50 mK, where 
qrTB ~ 1-96, the compensation is not strong, but the 
dephasing rate is already small here. The effect of the 
pressing field Ej_ remains smaller than from the polar- 
ization single-ripplon coupling for E± <C 10 3 V/cm. 
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Dephasing summary 

Ideally pure dephasing would be limited by the in- 
teraction of the electron motion with thermal ripplons. 
This would set a dephasing time of sa 0.5 ms even for 
T=100 mK. This value would be very sensitive to temper- 
ature and could be reduced significantly at temperatures 
lower than 50 mK, giving it a distinctive signature. It 
is possible, perhaps likely, that initial experiments will 
be dominated by anomalously high contributions from 
outside sources such as helium level fluctuations due to 
vibrations or spurious noise on the bias line, but given the 
analysis presented here it appears that it will be possible 
to eventually reduce these to acceptable levels. 

SPIN DYNAMICS 
Coupling to electromagnetic field 

The out-of-plane magnetic field B z that depends on 
the coordinates of in-plane motion (see Fig. 2) gives 
rise to spin-orbit coupling with Hamiltonian H s = 
—2/j,BS z xdB z /dx. Since the orbital motion is coupled 
to the cavity mode, with Hamiltonian H g = —exE [E 
is the operator of the electric field of the mode), the in- 
teraction H s b leads to a spin-mode coupling mediated by 
lateral electron vibrations. Formally, for weak coupling, 
this interaction is described by the spin-mode vertex cal- 
culated in one electron-vibration line approximation. It 
can be easily found also without diagrams. 

We will assume that the spin-orbit coupling is nonres- 
onant, so that \u x — ujl\ largely exceeds all relaxation 
rates and the detuning of the cavity mode from the Lar- 
mor frequency. A simple way to describe the spin-mode 
coupling is based on the equation of motion for operator 
x in the Heisenberg representation. If we disregard the 
spin-orbit coupling, this equation reads 

mi + mu? x x = eE. (27) 

We seek the solution that describes forced vibrations 
at the bare Larmor frequency ojl, which is very close 
to the cavity mode frequency. Then we obtain x « 
— eE/[m(uj\ — oj x )]. This solution is substituted into 
H s instead of x, which leads to the spin-mode cou- 
pling oc s z E. The coupling constant g s is given by 
Eq. (2) of the main text. We note that the spin-orbit 
coupling causes also renormalization of the Larmor fre- 
quency oc (fi B a x dB z /dx) 2 h' 2 uj L / (oj 2 l - wg). 



Spin decoherence 



In the absence of the enhanced spin-orbit interaction 
the lifetime of electron spin states is expected to ex- 
ceed seconds [2]. When the spin is coupled to the mo- 
tion, it will also inherit the orbital decoherence mecha- 
nisms studied above with a matrix element reduced by 
cx [iBd x B z a x /huj x . These mechanisms can be reduced 
by turning off the gradient field or changing the spin- 
motion detuning, to reduce the coupling. In addition 
to decoherence felt through the spin-orbit coupling, the 
electron spin can be dephased by fluctuating magnetic 
fields. One source of magnetic field noise is Johnson (cur- 
rent) noise in the trap electrodes, a white noise given by 
Si = 4fcT Aw/Re Z, where Aw is the bandwidth of in- 
terest, and the real part of the impedance is typically 
Rc Z <~ 50 SI. The field created at the electron by a 
current in the wire is B = 4 x 10 3 mG/mA, correspond- 
ing to a frequency shift of 5lul = 3.50 x 10 7 s _1 /mA. 
The Johnson noise is extremely small lOOpA/Hz 1 ^ 2 , giv- 
ing a coherence time X^j <~ 20 s. Another possible 
source of magnetic field noise is 1// flux noise [8]. of- 
ten seen in SQUID experiments with magnitudes of or- 
der S^,* = 10~ 6 $o/Hz 1/ ' 2 . The trap geometry con- 
tains no superconducting loops or Josephson junctions so 
this mechanism may not cause decoherence for an elec- 
tron on helium. If present it is reasonable to assume 
the flux would be distributed evenly over the trap area 
(~ 500 x 500 inn 2 ). Even in this conservative scenario the 
dephasing rate would be 1^$ <~ 200 s _1 , allowing many 
coherent operations. 
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